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Abstract We study the conformational dynamics within homo-polymer globules by solvent-implicit Brow- 
nian dynamics simulations. A strong dependence of the internal chain dynamics on the Lennard- Jones 
cohesion strength e and the globule size Ag is observed. We find two distinct dynamical regimes: a liquid- 
^s^j like regime (for e < e s ) with fast internal dynamics and a solid-like regime (for e > e s ) with slow internal 

dynamics. The cohesion strength e B of this freezing transition depends on Nq. Equilibrium simulations, 
where we investigate the diffusional chain dynamics within the globule, are compared with non-equilibrium 
C+H simulations, where we unfold the globule by pulling the chain ends with prescribed velocity (encompass- 

ing low enough velocities so that the linear-response, viscous regime is reached). From both simulation 
protocols we derive the internal viscosity within the globule. In the liquid-like regime the internal fric- 
tion increases continuously with e and scales extensive in Ng- This suggests an internal friction scenario 
where the entire chain (or an extensive fraction thereof) takes part in conformational reorganization of the 
globular structure. 
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1 Introduction 



Conformational dynamics of polymers play a crucial role in biological systems. For example, pore-translocation of 
I polymers such as RNA requires bond breaking and drastic conformational rearrangements to accommodate the geo- 
metrical constraints imposed by a pore [IH5] ■ Mechanical unfolding of biopolymers in force spectroscopy experiments 
induces profound conformational changes [4] > which typically involve dissipative effects [5j |6] . RNA sequences known as 
00 riboswitches experience conformational changes upon binding of small metabolites [7] . All of these transitions necessi- 
ty tate spatial rearrangements of the molecules and in certain cases require chain reptation within a collapsed region [B] . 
The time scale on which these changes happen is influenced by the medium, particularly the solvent viscosity, and by 
the polymer itself, through internal interactions. Therefore, if one desires to understand the dynamics of processes such 
as protein folding [5HT2"] , packing of DNA in the chromosome [TSHTo] , polymer collapse [T7HEJ ; or adsorption [2Tjjj2"2"] , 
t— I all dissipation and viscous effects have to be considered [25] . 

While the viscosity of the medium certainly influences polymer chain dynamics, there are a number of effects that 
together give rise to what is called internal friction. Local interactions such as conformational transitions of backbone 
bonds and dihedral angles [I2EI], entanglement effects and excluded volume interactions in polymer systems, degrees 
of freedom orthogonal to the reaction coordinate [25] , and the breakage and reformation of cohesive bonds [10l ESI EI] 
all lead to dissipation and thus increase the internal friction. For globular homopolymers, proteins in the molten 
globule phase [5S], and disordered intermediates during protein folding ES], these effects can be conceptualized as 
roughness on a hypothetical free energy landscape, corresponding to many competing and intermediate states [501 15T] . 
leading directly to the idea of an effective internal viscosity landscape [52Tl34| . There have been a large number of 
coarse-grained simulation studies on the force-induced unfolding of proteins [35] , globular polymers [17l ED I36H38] , 
and the diffusion of knots along a stretched chain [59fHI] . Cohesive interactions between polymer monomers have been 
shown to lead to a phase transition from a liquid- like to a solid-like globule for long enough chains [4"2Tf5U] . 

In this paper we study internal friction in two model systems, which both can be realized experimentally. We 
perform solvent-implicit Brownian dynamics simulations on a homopolymer. Attractive interactions are modeled with 
a Lennard- Jones potential, where the cohesive strength e is varied. First, we study the diffusion of a globule, which 
forms from a polymer chain held at constant extension smaller than the contour length. This simulation is conducted 
under equilibrium conditions and no external forces are applied. For e < e s , where e s depends on the globule size, 
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pronounced diffusion of the globule relative to the linker chain section is observed and characterized by the globule 
diffusivity Dq. Dq is a direct measure of the internal globule viscosity, since motion of the globule relative to the 
linker chain requires internal rearrangements. We observe that the diffusion constant is proportional to the reciprocal 
globule size Dq ~ Nq 1 and shows a marked dependence on the cohesive strength e. For e > e s , no diffusion is observed 
on the simulation time scales and the globule is stuck in a single conformation. This reflects the change of the internal 
dynamics going from liquid-like (e < e s ) to solid-like (e > e s ). Support for this interpretation is given by our second set 
of simulations, where we measure the dissipated work during stretching and unraveling the globule with a prescribed 
finite velocity. This is the same setup as studied by us before, but here with significantly longer chains and slower 
pulling velocities |32j . Like in our equilibrium simulations, we observe two different regimes. For large e > e s , the 
force extension curves are characterized by pronounced fluctuations, which are absent for e < e s . Pulling decreases 
the number of monomers inside the globule, causing the fluctuations to vanish once Nq is below a certain threshold. 
Therefore, reducing the number of monomers Nq in the globule by stretching the chain, drives the system from the 
solid into the liquid state, similar to the reduction of the glass temperature of polymers close to interfaces [5TH54"] . In 
the liquid state we perform extensive simulations and show that the internal viscosity is extensive, meaning that the 
dissipated work per monomer that is pulled out from the globule scales linearly with the globular monomer number Nq 
and the pulling velocity v, similar to our findings from the equilibrium simulations. For the chain dynamics this means 
that in the liquid-like regime, a model based on local viscous friction is valid, but that the viscosity depends on the 
size of a globule. Clearly, internal friction effects dramatically influences the time scales of chain dynamics in globules. 
In the solid-like regime, no definite conclusion on the conformational dynamics is possible from the simulations we 
performed. 



2 Brownian dynamics simulations 

We model the homopolymer by N freely jointed beads of radius a, interacting with a potential [/({r^}), which 
depends on the set of positions of all beads {rW}. The position of the i th bead obeys the overdamped Langevin 
equation [551 I5S] 

— = -rt)V P «,tf({rW}) + Mo/ W (t) ) (1) 

where (iq — l/(67r?7oa) is the Stokes mobility of a sphere with radius a in a solvent with viscosity i]q. is the 

random force acting on the i th bead. The components of the random force /£ (t), a = x, y, z, satisfy the Stokes- 
Einstein relation 

(^ ) W/? ) (0) = — Sii&ctW-t) (2) 
Mo 

and are unbiased {fa if)) = 0. Hydrodynamic interactions are neglected, since we are interested in the internal friction 
of the polymer caused by the monomer interactions. We express energies in units of the thermal energy keT, lengths 
in units of the bead radius a, and times in units of t = a 2 /(/iokeT), which is the characteristic diffusion time of a 
single bead. Using the dimensionless quantities t — t/r, f — r/a, f = f /{kg,T /a), and U — U/{k^,T), eq. ([I]) reads 

° w = -V fW U{{f^}) + p\t). (3) 

For the Brownian dynamics simulation we use a discretized version of eq. ^ with a time step At = 0.0005r and 
obtain the coordinates at time t = nAt [57] 

fW(n + l) = fW(n) + (-V.ot/flfW}) + / w (n)) At (4) 
and 

(f a i) (n)f^(n')) = ^S ij 5 a0 S nn> . (5) 



3 Diffusivity of a globule along a periodic chain: equilibrium simulations 



In our preceding work we determined the internal viscosity t/q of a homopolymeric globule by measuring the dissipated 
energy when unfolding the globule by pulling apart the chain ends at finite speed [32) . We showed that for moderate 
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Figure 1: The diffusion of a globule relative to the stretched polymer linker is simulated. The globule mobility, the 
number of monomers Nq inside the globule, the diameter of the globule Iq, and the fluctuations of the linkers (i.e. the 
stretched part of the chain that is not part of the globule) depend on the cohesive strength e, eq. (10 1. Snapshots at 
different times for e = 0.8 (left) and 2.09 (right) are shown. To prevent motion of the linker chain section and thereby 
obtain directly the globule diffusivity relative to the chain, the monomer i tr in the middle of the linker chain (indicated 
by the arrow) is trapped by a harmonic potential, eq. (11 1. Periodic boundary conditions for a box of length L using 



the minimum image convention are employed to model an infinite polymer. The color coding indicates the running 
monomer index along the chain. 



cohesion e the internal friction is the dominant dissipative effect as long as the majority of the monomers are part 
of the globule. However, as the globule is unraveled, the size of the globule decreases and more and more beads in 
the linker sections dissipate energy due to solvent friction. These effects had to be subtracted in order to obtain the 
internal globule friction. Likewise, reaching the relevant linear-response regime at low pulling velocities is subtle. Here, 
we introduce a novel system setup, where we perform equilibrium diffusion simulations with a globule of constant size: 
The entire polymer chain is held at constant extension smaller than the contour length, so that a globule forms for 
high enough cohesive energy. Since we fix the position of the linker chain in space, as we explain further below, the 
diffusive motion of the globule necessitates internal rearrangements, so that the globular diffusivity relative to the linker 
chain is a measure of the globule internal friction. Since the method works at equilibrium, the linear viscous regime 
is automatically obtained. By this method, the internal viscosity inside the globule manifests itself as a macroscopic 
and experimentally observable quantity: the diffusivity of the globule, Dq ~ 1/?7g- This system could be realized 
experimentally by a polymer held at constant extension in an optical or magnetic tweezers setup |58j . In such an 
experiment the cohesive force can be varied by changing the solution conditions, whereas the size of the globule can 
be varied by changing the trap distance. 



3.1 Model 

3.1.1 Description of the system 

We consider a polymer held at a fixed extension L = 50a in the x-direction, which is smaller than the contour length 
L c = 2a (N — 1) of the polymer. A globule will form for large enough attractive Lennard- Jones interaction between the 
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monomers, see fig. [T] To eliminate finite size effects, we introduce periodic boundary conditions in ir-direction, which 
are implemented via the minimum image convention |59j : The components of the vector pointing from r'*' to are 
given by 

= 2?(r« , r« ; L) = ((»#> - r« + 3L/2) mod L) - L/2 (6a) 
(y) = Cj) _ W (6b) 

We use a box size L = 50a in all simulations. The potential energy has four contributions 

U=U h + U r + U LJ + U tT . (7) 

Uh and U-[ are the bond potentials acting between neighboring monomers. The backbone bonds are modeled by 
harmonic potentials 

N-l 
i=l 

with k = 200kBT/a 2 and r^' J ' — |r' l ' J ')|, see eq. ([6]). As periodic boundary conditions are employed, the polymer forms 
a closed ring, which is achieved by connecting the first and last monomer by 

U I = ^r^-2ar. (9) 



The monomer cohesion and excluded volume interactions are modeled with a Lennard- Jones potential 

^-'sgtur-'U)*)- 

e = models an ideal phantom chain without excluded volume interactions and attractive interactions between the 



monomers. For e > the first term in cq. ( f ) accounts for the repulsive excluded volume interaction at short separa- 
tions, whereas the second term is responsible for cohesion which reflects hydrophobic attraction between monomers in 
a solvent-implicit fashion. For < e < £q ~ 0.5 the polymer is in the swollen state and no globule exists. Increasing e 
above eq causes the polymer to collapse and a globule forms. The static globule behavior from our simulations agrees 
with previous work [T71 l3"2l r4"2"H5"0l I^U] . The value of the cohesive strength at the globule transition, eq, depends on 
the system size [43]. For even larger values e > e s a solid phase appears. 

As we are interested in the relative motion of the globule with respect to the rest of the chain, we prevent the 
linker chain from moving by an external trapping potential acting on one monomer, £/t r , which can be viewed as the 
effect of an e. g. optical tweezers. The linker chain is the stretched part that does not belong to the globule, fig. [T] We 
introduce a harmonic trap potential 

U tr = ^(r^-Rr, (11) 

which is located at R and acts on bead iti that is in the middle of the linker, with stiffness n tr — lOkeT '/a 2 . The index 
of the trapped bead i tr depends on the position of the globule and on the set of beads that belong to the globule. 
The exact definition of the globule is described in the next section. Due to globule motion, the index of the bead in 
the middle of the linker might change from i tr at time t to i[ T at some later time t' . If this happens, we update the 
^-coordinate of the trap position to the new position R x (f) = R x {t) + 2aT>(i tT , i' tI ; N), eq. and the trapping force 
then acts on bead i' tI . The trapped beads in the snapshots shown in fig. [I] are indicated by arrows. 

This setup allows us to study the diffusion of the globule relative to the linker chain along the x-axis. The effect 
of using periodic boundary conditions rather than long linkers to each side is fourfold. First, the simulation is sped 
up as the system is smaller. Second, the globule is stabilized as the configurational space for the unfolded system is 
reduced since only fluctuations up to a wavelength of the order of the simulation box are possible. Third, motion of 
the globule as a whole without internal friction is reduced. Fourth, knots may not form as the polymer forms a closed 
ring. 



3.1.2 Definition of the globule 



We partition the a:- axis into bins of width = 2a and measure the monomer density pk in the fc th bin, k = 1, . . . , L/l\,. 
The globule is defined as the region, where the monomer density projected on the cc-axis fulfills the condition pk > 5/a. 
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Figure 2: The density profile of the globule depends on the cohesive strength e. The globule is defined as the region 
where the monomer density is p > 5/et (indicated by the horizontal line). 



In fig. [2] the density profiles of the snapshots in fig. [T] are shown. If at the edge of the globule the density profile is not 
monotonous and pk < 5/a, but Pk-i-.Pk+i > 5/a, we account for such cases by adding these bins to the globule, too. 
This ensures that we end up with a list of bins that are connected. For small e < 1 and small N < 150 it is possible 
that - according to the above definition - more than one globule exists or that the globule is smeared out over the 
complete simulation box, however we will not consider simulations in which this occurs. Such complications are never 
observed for e > 1 and TV > 150. 



3.1.3 Definition of the index of the trapped bead 



Let us first define the index of the bead i T , which is at the right edge of the globule, and the index of the bead i\, 
which is at the left edge of the globule. i r is obtained by picking a monomer inside the globule and moving along the 
chain contour with increasing monomer index. i r is the largest index that is still in a bin belonging to the globule. 
We also check for loops, which leave the globule and return again: If a loop occurs, we add all monomers of the loop 
to the globule even if they lie in a bin outside the globule. i\ is defined in the same way yet by decreasing the index. 
The index of the central bead i c in the middle of the globule and the number of monomers Nq inside the globule are 
defined as i c = (i\ + i r )/2, Nq = i T — i\ + 1 if i\ < i r and i c = {i\ — N + i r )/2, Nq = i\ — i T + 1 if i\ > i r . The index of 
the central bead i c yields the index of the trapped bead via i tI = i c + N/2, which therefore depends on the motion of 
the globule. We could have also trapped a fixed monomer for the whole duration of the simulation. Our procedure (i.e. 
trapping a bead that by construction is never part of the globule) allows maximal dynamic freedom for the globule 
and thus improves equilibration of the system. 



3.2 Results 



3.2.1 Number of monomers inside the globule 



The number of monomers Nq inside the globule depends on the chain length N and on the cohesive strength e, fig. [3^,. 
Increasing e raises Nq only up to a limiting value Nq, which is determined by the equation 

(iVa) 1/3 + (N - N G ) = L/(2a) . (12) 



Eq. ( |12| describes a close-packed spherical globule consisting of N G monomers with a tightly stretched linker consisting 
of N — N G monomers. As the limiting value of our simulations coincides with the predicted value of eq. ( 12 1, as shown 



in fig. [3^, our definition of the globule is justified. For large e, the linkers are completely stretched, as the energetic 
gain of a monomer joining the globule outweighs the entropic loss of reducing the fluctuations of the linker. For small 
e, monomers are not tightly bound to the globule as can be seen qualitatively from fig. [T] Therefore, size fluctuations 
of the globule are more substantial for small cohesive strengths and decrease upon increasing e, as shown in fig. [3b. 
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Figure 3: (a) The number of monomers Nq inside the globule increases weakly with the interaction strength e until 
it finally levels off. The limiting values N G are depicted by horizontal broken lines and are determined by eq. (12), 
describing a spherical globule and tightly stretched linkers. Error bars denote the standard deviations of the distribu- 
tions, (b) Probability distribution of the number of monomers Nq inside the globule for chain length N — 200 and 
different cohesive strengths e = 0.8, 1, 1.5, 2.81. For small e considerable size fluctuations are observed, which decrease 
upon increasing s. Further, increasing e raises the mean Nq up to the limiting value Nq given by eq. (12) (broken 
vertical line). 



3.2.2 Definition of the center of the globule 

In order to properly describe the motion of the polymer in a simulation box with periodic boundary conditions in 
x-direction, we introduce non-periodic polymer coordinates {rW}, which are not restricted to the primary simulation 
cell. They are defined recursively starting from the first monomer with = via 

^(i+i) _ + r (M+D 5 ( 13 ) 

for i > 1, where r^' (+1 ^ is given by eq. ([6]). Since we are considering a ring-like polymer with periodic boundary 
conditions in the x-direction, the center of mass of the globule has to be defined carefully. The center of the globule 



is calculated by using the non-periodic polymer coordinates, eq. (13), according to 

-(G) - 



, = — Y>W . (14) 

Nq^ 1 ' 



Periodic boundary conditions may introduce jumps of the size of the box L in the trajectory r^(t). We remove 
those jumps by connecting the value of any quantity q(t) at time t, which is subject to periodic boundary conditions, 
with the value q(t — At) in the previous time step via the minimum image condition. 



3.2.3 Diffusivity of the globule 

Our setup enables us to study the motion of the globule in space in a fashion that is coupled to its internal confor- 
mational dynamics. In other word, since we fix the linker chain in space, only internal chain rearrangements within 
the globule lead to diffusion of the globule in space. In fig. H trajectories of the x-coordinate of the center of the 
globule, eq. (14), are shown for N = 200 and various cohesive strengths e. Fig. H demonstrates that the diffusivity 



decreases with increasing e. As the monomers become more cohesive, it is more dTmcult for the globule to rearrange 
internally and hence to move. This effect is further accented by the increase in globule size Nq, see fig. [3J which 
additionally decreases the mobility of the globule. 

To quantify these observations, we calculate the mean squared displacement (MSD) of ri G ' . For normal diffusive 
behavior one expects the MSD to scale linearly with time and to be characterized by the diffusion constant Dq 

MSD(t) = ((rf\t) - ri G )(0)) 2 ) = I f {r^\t + f ) - r^\t')) 2 dt' = 2D G t . (15) 
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Figure 4: Trajectories of the x-coordinate of the center of the globule, r x , defined in eq. (14), for N = 200 and 
various interaction strengths e — 1. 1.25, 2.5, 2.81. The smaller the attractive interaction between the monomers, the 
more mobile the globule is. If e is large, the globule is frozen in a single conformation and does not move at all. 




Time t j i 



Figure 5: The mean squared displacement (MSD) of the center of the globule for N = 200 and e — 1, 1.75, 2.5, 2.81 
is calculated from the trajectories, fig. El using eq. (15). Symbols denote the measured MSD from our simulations, 
lines the corresponding linear fits for MSD(i) > 10. For e < e s with e s w 2.3 for N — 200, the globule exhibits normal 
diffusion and the diffusivity Dq decreases as e increases. For e = 2.5, normal diffusive behavior is obtained, whereas 
for e = 2.81 the diffusion time scale is of the order of the simulation time and normal diffusive behavior is barely 
reached. 



MSD curves for N — 200 and various cohesive strengths are shown in fig. [5] on a double logarithmic plot. The MSD 
curves are fitted with linear functions for MSD(t) > 10 in order to obtain the diffusion constant Dq. For small e < e s , 
normal diffusion is observed with Dq decreasing as e increases. However, as can be seen in fig. [5] for TV = 200 and 
e = 2.5, the normal diffusive regime only occurs at very long time scales and for e = 2.81 is barely reached on the time 
scales of our simulations. We attribute this to a change of the internal dynamics of the globule, which is slowed down 
with increasing cohesion and is effectively suppressed for e > e s . One remaining pathway for the globule to rearrange 
at very large cohesion is to dissolve - at least partly - and refold into a different configuration. As a consequence, 
the time scale characterizing the internal dynamics should become comparable to the time scale on which the globule 
dissolves. This dissolution time is huge as it scales exponentially with eNq and is beyond our simulation time. For that 
reason we observe only stuck globules for large cohesive strengths, which remain in a single conformation. The fitted 
diffusion constants are shown in fig. [6] and contrasted with an idealized limit, where internal friction is absent and the 
globule and the linker move independently. The diffusivity in this limit is given by the Rouse diffusion constant of N G 
monomers 1551 



(16) 
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Figure 6: (a) Diffusion constant Dq of the globule (open symbols) as obtained from linear fits to the MSD curves, fig. [5] 
Dq decreases with increasing N and e. At e = e$, with e s depending on N, a transition of the internal dynamics occurs 
and the diffusion constant drops to zero. Dq is compared to the ideal Rouse diffusion constant Dq (solid symbols) of 
a globule with N' G monomers, eqs. (16) and (17), which can move freely and independent of the linkers, (b) Rescaling 



the actual diffusion constant by the ideal diffusion constant, Dq/Dq, removes the N dependence to a large extent in 
the liquid regime. 



N'q < Nq is the reduced number of monomers within the globule that actually have to comove when the globule is 
displaced by some distance x. N' G can be estimated by continuing the linker through the globule and subtracting the 
number of monomers, which belong to this internal linker section, from Nq, see fig. [7] for an illustration, 

N g = Nq-(N -Nq)— ^- . (17) 

L — Lq 

In fig. [6^, a pronounced dependence of the diffusion constant Dq (open symbols) on the globule size Nq is observed. 
The diffusivity decreases as the internal interactions e increase. This is due to the coupling of the internal dynamics of 
the globule to the overall motion of the globule, because, as explained before, in our simulation setup the linker is fixed 
and the globule can only move via internal conformational chain reorganization. Slowing down the internal dynamics 
by increasing e thus reduces the mobility of the globule. For small e one observes Dq sa D (filled symbols) , implying 
that internal friction is unimportant. Conversely, increasing e causes increasing deviations between the mobility of 
the globule and the ideal system, until finally Dq drops to zero as the liquid-solid transition is crossed. For e > e s 
the globule is in a frozen state with strongly suppressed internal dynamics. Since the linkers are trapped, this also 
impedes the motion of the globule as a whole and leads to a vanishing diffusion constant. To disentangle size from 
cohesive effects, we show the rescaled diffusivity Dq/Dq in fig. ^jp. For e — > the rescaled diffusivity approaches 
unity, indicating that internal friction is unimportant. The rescaleddiffusivity exhibits for small cohesion only a weak 
dependence on the size of the globule. As s approaches the solid regime, however, deviations between different system 
sizes become observable. Thus, we propose that in the liquid regime, the internal friction is extensive and to leading 
order scales with the size of the globule 

D G ~ l/N' G (18) 

as will be corroborated by our non-equilibrium simulations in section [4] Further below we will also show how the 
internal viscosity within the globule can be extracted from the diffusivity ratio Dq/D . 

4 Forced unfolding of globules: non-equilibrium simulations 
4.1 Model 

We now discuss our simulations for the non-equilibrium force-induced dissolution of homopolymeric globules. Here, 
the two chain ends are positionally constraint at r' x '(t) = — R(t) and r^ N \t) = R(t). The trap positions are moved 
at constant speed R x (t) = i?i mm ' + vt with R y / Z — 0, from iii m ^ up to a maximal position i?i max ' 1 = a(N — 1), 
producing a time- varying chain extension x(t) = 2R x (t). During the pulling, the force acting on the terminal beads 
is measured. The backbone bonds are modeled by a harmonic potential Ub as before, see eq. ([8]). Excluded volume 
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Figure 7: Illustration of the concept and calculation of the reduced number of monomers N G in the globule. For a 
displacement of the globule by a distance x only a reduced number of monomers have to move along, here N G = 14 
(shown in orange) . In order to obtain N G , the linker is continued thr ough the globule and the monomers belonging to 
this internal linker are subtracted from N& yielding N G , see eq. (17). 



and cohesive interactions are again modeled by a Lennard- Jones potential JTlj, eq. (10). Eq. Q is used to integrate 
the Langevin equation. We no longer employ periodic boundary conditions, however still prevent knot formation by 
introducing a potential, which mimics two repulsive bars that extend from the first/last bead to the left /right along 
the rc-axis 



2 (2a/ p 



(i)) 



lj if pW < 2a and \ri l '\ > R 
else, 



(19) 



Sir 



p\-i = y r y ' +ri' ■ does not affect the stretching response. The total energy is therefore given by U = 
Ub + f^Lj + U^. Two different protocols, annealed and un-annealed, are used for the initial configurations in order 
to investigate the history dependence on the globule pulling response. We record force extension curves for vari- 
ous chain lengths N = 50, 100, 200, 300, cohesive strengths < e < 4.1, and pulling velocities v — v/(a/r) = 
0.001, 0.0045, 0.01, 0.0225, 0.045 going significantly beyond our previous work, where the largest system was N = 100, 
the slowest velocity was v = 0.0045, and only un-annealed initial configurations were used |32j . For each parameter 
set, twenty stretching cycles are simulated. 



4.1.1 Annealed initial structures 

Annealed initial structures are obtained by performing an equilibrium annealing simulation with moderate cohesive 
strength e = 0.8 and fixed trap position R x = O.laN. The LJ interaction is strong enough to induce globule formation, 
yet small enough to allow for rapid equilibration of the chain conformation inside the weakly collapsed globule. Knot 
formation is prevented by virtue of the potential in eq. (19). Every t = 20000r a structure is recorded, which is 
subsequently equilibrated for t — 10000t using the target cohesive strength at which the pulling simulation is to be 
conducted. The resulting structure is used as one initial configuration for the subsequent pulling cycle. In fig. [8^, the 
pulling protocol is illustrated and in fig. [9] typical initial configurations are depicted. 



4.1.2 Un-annealed initial structures 

This set of initial structures is obtained by starting from an extended configuration and moving the traps from 
i?i max) = a(N — 1) to R x = QMaN. Typically, we observe the formation of one globule for v < 0.01 and the formation 
of two globules near the traps for v > 0.01, which merge at small extension, see snapshots in fig. [8]b. Without pausing, 
the traps are subsequently extended to R x = O.laN , where the actual pulling cycle starts and force extension curves 
are recorded for further analysis. During the compression stage, the traps are moved with the same velocity v with 
which the force extension curve is recorded, fig.[8]D. For larger cohesive strengths e, non-equilibrium ordered structures 
prevail as initial configurations for the subsequent pulling cycle, see fig. [9] for an illustration. Although not the prime 
target of our present investigation, such structures might be of relevance when studying the dynamics and packing of 
DNA chromatin structures far from equilibrium [13H15j . 
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Figure 8: Illustration of the two different pulling protocols. R x and — R x denote the positions of the two traps that 
positionally constrain the chain ends, (a) Preparation of the annealed structures. After a long equilibrium run with 
e = 0.8 and R x =0.1 (broken line), the globule is equilibrated at the cohesive strength at which the pulling curve is 
recorded (horizontal solid line). The subsequent pulling cycle is depicted by the thick gray line, (b) Preparation of the 
un-annealed structures. The traps are moved from complete extension, R x = a(N — 1), to R x = 0.03aN, which ensures 
that one single globule forms. Without pausing, the pulling cycle starts and force extension curves are recorded and 
analyzed in the interval 0.1 < R x /(aN) < 1 (thick gray line). 




Figure 9: Typical initial configurations in the annealed (top) and un-annealed simulations (bottom) for e = 2.91 and 
v = 0.001. For the large e shown here, the collapsed un-annealed structures show some residual ordering. The color 
coding indicates the monomer index along the chain contour. 



4.2 Liquid-solid transition for large cohesive strengths 

In fig. |10| stretching curves averaged over twenty pulling cycles for the annealed pulling protocol are shown for various 
cohesive strengths e and v = 0.001, iV = 300. Beyond the globule transition, e > eq 0.5 [32J HH [SU] , an £ dependent 
force plateau is observed. The plateau force F p increases as e increases and is, for the relatively slow pulling speed 
shown here, mostly associated with the equilibrium free energy per unit length of globule formation, F p a ~ e — sq |32j . 
For large extensions of the order of x/L c « 0.8, a dip in the force extension curve appears, which is the signature 
of a pulling induced globule dissolution. For even larger extension, essentially no monomer contacts are present, and 
the force extension curve becomes independent of the cohesive strength e and follows the trace of an extensible freely 
jointed chain (shown as a broken line) 

x/L c = coth(2aF/(k B T)) + k B T/(2aF) + F/{2an) . (20) 



A phantom chain (e = 0) coincides perfectly with eq. ( |20| ). The curve for e = 2.5 features a maximum at x/L c w 0.3, 
whose origin will be discussed in the next paragraph. 

Fig. [IT] shows force extension traces for the two different pulling protocols for N — 300, v = 0.001, and relatively 
strong cohesive forces e = 2.08, 2.5, 2.91. The thick curve in each plot depicts the average over all 20 pulling curves, 
while the thin curves show individual force-extension traces for different initial configurations. In figs. |ll^ -c pulling 
curves starting from un-annealed configurations are shown. Increasing e leads to an increasing plateau force similar 



to fig. 10 For e > 2.5 stronger fluctuations of the force are observed, but the force extension traces are qualitatively 
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Figure 10: Force extension curves for constant velocity v = 0.001, various cohesive strengths and ./V = 300 using 
annealed initial configurations. All curves are averages over twenty pulling cycles. Above the globule transition, e > 
Eg ~ 0.5, a force plateau followed by a dip in the force extension curve is observed. The broken line depicts the 
theoretically expected force extension trace of an extensible freely jointed chain for e = 0, eq. (20). For e = 2.5 the 
pulling curve exhibits a marked maximum at small extensions, which indicates frozen internal dynamics. 



similar to the curves for s = 2.08. The situation is vastly different for pulling curves starting from the annealed initial 
configurations, figs. [TTpl-f. Again, for small cohesive strengths e < 2.08 the pulling curves are smooth and no strong 
fluctuations occur, cf. fig. |10| However, increasing the cohesive strength further leads to pronounced fluctuations in the 
force extension curves. This is due to a transition of the internal dynamics from liquid-like to solid-like [3"2l l4"2l B51 Pf81 - 
[50) [61] . Since the un-annealed initial configurations are rather ordered - especially for large £ - the globules are easily 
unwound by simply retracing the configurational intermediates that were encountered upon folding in reversed order. 
The liquid-solid transition occurs for the un-annealed structures, too, but has almost no effect on the non-equilibrium 
pulling simulations. This is very different for the annealed simulation protocol, where the mean force is higher and the 
variation of individual force curves around the mean force is also pronounced. 

The dependence of the liquid-solid transition on the globule size Nq is illustrated in fig. [12J Here, force curves 
obtained using the annealed pulling protocol are plotted versus the number of monomers Nq inside the globule for 
different monomer numbers N and cohesive strengths. Nq is calculated via a modified version of eq. (12) where the 
extension x of the polymer is used instead of L. One notices that the pronounced noise in the curves ceases once the 
globule is below a certain size, i.e. once Nq is below a certain threshold value, where this critical size depends on e. 
This feature is independent of the chain length N but solely depends on Nq, i.e. curves with equal e but different N 
coincide once the globule size has fallen under the critical size |4"2Tf5"0"] . The cohesive strength e s at which this liquid- 
solid transition occurs is also independent of the pulling velocity within the range of velocities studied. In fig. 13 we see 



that for e = 2.08 even with the highest pulling velocity, no large fluctuations in the force extension curves are induced 
and the globule remains in the liquid phase. Therefore, the liquid-solid transition at e s is not a mere non-equilibrium 
pulling feature but an indication of a change in the equilibrium internal dynamics. The transition observed in the 
non-equilibrium pulling simulations is of course the same transition that induces the abrupt change of the equilibrium 
globule diffusivity, demonstrated in fig. [6j since both the equilibrium globule diffusion and the non-equilibrium pulling 
response are ultimately related to the configurational chain dynamics inside the globule. 



4.3 Internal friction 



4.3.1 Internal friction is viscous 



For small e < e s , the internal chain dynamics is fast and no significant difference between the simulations with annealed 
and un-annealed initial configurations is observed. In fig. 13 pulling curves for various pulling velocities v are shown. 
We observe that even for a cohesive strength as large as e = 2.08 and polymers as long as N = 300, the two slowest 
pulling curves v = 0.001, 0.0045 almost coincide. This indicates that force extension curves for v = 0.001 are already 
very good approximations for equilibrium curves. Increasing velocity leads to increasing energy dissipation. There 
are two major mechanisms leading to dissipation: solvent friction and internal friction. These dissipation mechanisms 
dominate in different parts of the pulling curve. For small extensions, when most of the monomers are part of the 
globule, internal friction dominates. Towards the end of the pulling process, the globule is markedly smaller and the 
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Figure 11: Pulling curves for N = 300 and v — 0.001 with (a-c) un-annealed and (d-f) annealed structures as initial 
configurations. The thin lines are individual force extension traces, whereas the thick black line is the average over 
twenty pulling curves. The un-annealed structures exhibit rather smooth pulling curves and no drastic differences 
between the various cohesive strengths e is observed. In contrast, the annealed structures feature marked fluctuations 
above a certain threshold e s ss 2.1 for N — 300, which diminish for large extension, i. e. when the globule sizes decrease. 
We attribute these fluctuations to a liquid-solid transition, which dramatically changes the internal dynamics (see text). 
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Figure 12: Averaged force extension curves in the vicinity of the liquid-solid transition are shown as a function of 
the number of monomers inside the globule, Nq, eq. (12 1, for the annealed set of initial configurations. Once Nq is 
below a certain threshold (« 100 for e — 2.91, ~ 140 for e — 2.5) the globule is driven into the liquid state, and force 
curves for different N and equal e collapse. The dependence of the liquid-solid transition on Nq is discussed in detail 
in ref. 021. 



solvent friction, which acts mostly on the linker chain, dominates the pulling curve. The dissipated work AW(v, e) is 
defined as the difference between the total work at finite velocity and the equilibrium work [3"2"] 



AW(v,e) =W(v, e)-W cq (e) , 
where the work done by one trap is generally defined as 



(21) 



W 



see the inset of fig. 
W(v, e) to v -t 0. 



F{x)dx , (22) 
for an illustration of this definition. The equilibrium work W cq (e) is obtained from extrapolating 



fig. 14 we show the dissipated work per monomer, AW/N , as a function of the velocity v for 
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Figure 13: Averaged pulling curves for various pulling velocities with N = 300, e — 2.08, and annealed initial config- 
urations. The friction force increases with increasing pulling velocity. The curves with the slowest pulling velocities 
v = 0.001, 0.0045 almost coincide, indicating that with these slow velocities one has approximately reached the equi- 
librium pulling limit. The dotted curve is one individual pulling curve for v = 0.045. The inset illustrates the definition 
of the dissipated work, eq. (21), which is the shaded area between the two curves. 
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Table 1: The critical cohesion strength e s of the liquid-solid transition depends on the length of the polymer N. For 
s < e s , the system is in the liquid phase and fluctuations of the dissipated work are small. For e > e s pronounced 
fluctuations are observed due to slow internal dynamics, see fig. 



1G 



N = 300 and various e. Below the liquid-solid transition, e < e a « 2.1 for N = 300, the dissipated work scales 
linearly with v, which is further demonstrated in fig. 15 where the ratio AW/ (vN) is plotted. This shows that in 



the liquid state, the simulations are conducted in the experimentally relevant linear response regime and the friction 
is essentially of viscous nature. Above the liquid-solid transition this scaling breaks down, or, in other words, the 
velocities probed in the simulations are not low enough to reach the viscous regime. Further, we observe that below 
the globule transition, e < e G ps 0.5, the dissipated work is almost independent of the cohesive strength. This suggests 
that monomer-monomer attraction and topological constraints, e. g. entanglements, are negligible for small values 
of the cohesive strength, i.e. in the non-collapsed state. In order to extract such subtle friction effects, much longer 
simulations would be needed. Knots, which in principle might arise due to bond crossing, are not observed in our 
simulations. 

The fluctuations of the dissipated work can be conveniently used to define the liquid-solid transition. In fig. |16| the 
standard deviation 



= y/k/(k-l) y/{AW 2 - {AW) 2 ) 



(23) 



obtained from k = 20 measurements of AW for v = 0.001 in the annealed pulling protocol is plotted against e on 
both linear and logarithmic scales. For e < e s , a aw /N is independent of N and e. As one crosses from the liquid into 
the solid phase at e — e s , the fluctuations increase. Fig. |l"6| clearly demonstrates that e s depends on the system size 
N. The estimates of e s extracted from the behavior of oawIN are shown in tab.JT] 



4.3.2 Internal friction is extensive 

We demonstrated that in the liquid state the dissipated work AW and thus the internal friction scale linearly with 
the velocity v. Therefore, we can express the friction force ff r in terms of a velocity independent viscous friction 
coefficient r(e, Nq), which depends on the number of monomers in the globule, 



F b (v,e,N G ) = F(v,e,N G ) - F cq (e,N G ) = v r(e,N G ) 



(24) 



where F eq is the equilibrium force extension curve. Analogous to the Stokes friction of a sphere, we define the friction 
coefficient as 



r(e,N G )^ar, G (e)N^, 



(25) 
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Figure 14: Dissipated work per monomer AW/N as a function of the pulling velocity v for N — 300 and different 
cohesive strengths e starting from the (a) annealed and the (b) un-annealed set of initial configurations. Symbols depict 
simulation data, lines show linear fits to the data according to eq. (27). Below the globule transition, e < eq « 0.5, the 
curves collapse, i.e. interaction contributions to friction are negligible. Below the liquid-solid transition, e < s s ~ 2.1 
for N = 300, the data is linear in v in the whole range of velocities studied, indicating that friction is of viscous nature. 
Above the liquid-solid transition the linear scaling breaks down. AW is slightly lower for the simulations starting from 
the un-annealed set of initial configurations showing that the globule is not completely equilibrated and still rather 
ordered. 



p 300 



S 200 


M 
O 

s 

■a 

s ioo 



(a) 



_ 1 1 1 1 + 1 1 1 I 


| I 1 

N 


i i : i : i 
= 300, annealed 




H — 1 

□ e 

— 0s 

- x — ° e ~ 


= 
= 0.4 
= 0.6 




- V e = 1.25 

O e = 2.08 

X e = 2.5 

+ e = 4.1 




X 

Ur-i 


X 




+ : 

X 




— — 




-o— 




[ 1 ! 1 1 1 ! 1 ! 1 1 ! 1 ! 1 ! 1 ! 1 1 1 ! 



0.01 



0.02 



0.03 



0.04 



Velocity v / (o/t) 



< 

■a 

o 



0J 
4-1 
CO 



400 
300 
200 
100 



(b) 



JV = 300, un-annealed 



□ e = 
■Oe = 0.8 
■ Ve = 1.25 
+ 



-Oe = 2.08 
-X e = 2.5 
- + £ = 4.1 



+ 

X 



0.01 0.02 0.03 

Velocity v / (o/t) 



0.04 



Figure 15: Dissipated work per monomer rescaled by the pulling velocity, AW/(Nv), as a function of the pulling 
velocity v for N — 300 and different cohesive strengths e starting from the (a) annealed and the (b) un-annealed set 
of initial configurations. Symbols denote simulation results, lines depict _T as obtained from linear fits according to 
eq. (27) in fig. 14 For e < e s , our simulations are carried out in the linear viscous regime and there is no v dependence 
after rescaling\ — Sbove the liquid-solid transition the linear scaling in v breaks down. For e = and 2.08, error bars 
indicate the standard deviations of the twenty measurements of AW(v,s). 



where J7g( £ ) is the internal viscosity, which depends on e but not on Nq. The exponent 7 describes the dependence 
on the monomer number inside the globule, and characterizes the friction mechanism at work during unraveling the 
globule. Two limits can be distinguished: first, if 7 = 0, the friction force is independent of the globule size Nq, and 
only a finite number of monomers, which does not scale with Nq, contribute to dissipation. We call this limit local or 
intensive friction, and the pictorial mechanism for such a scenario could be that monomers are peeled off one-by-one 
from the surface of the globule. Second, 7 = 1 describes the limiting situation where a finite fraction of the globule 
that is proportional to Nq, or even the entire chain, rearranges in the unfolding process and hence contributes to the 
friction force. We call this limit global or extensive friction and note that there are several conceivable microscopic 
mechanisms. One possibility could be a reptation scenario, where a reorganization of the globule conformation is 
accompanied by a reptation of the entire globular chain section (or a chain fraction of length proportional to Nq) 
through the globule. Another possibility would be that the stress propagates through the whole globule and leads to 
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Figure 16: The standard deviation per monomer (Jaw/N of the dissipated work, as obtained from twenty pulling cycles. 
For e < e s , (Jaw /N is independent of the chain length N and small. The fluctuations of the dissipated work increase 
significantly when e > e s , where e s decreases for increasing JV, see tab. [I] The data is obtained from simulations with 
v = 0.001 and annealed initial configurations. 




■ON = 50 
• □ N = 100 
■ N = 200 
-A JV = 300 

— e = 

- e = 1.25 

— e = 2.08 

- e = 2.5 

— e = 4.1 



0.04 



Velocity v / (a/f) 



Figure 17: Rescaled dissipated work per monomer AW / (N 2 v) as a function of the pulling velocity v for different JV 
and cohesive strengths e starting from the annealed set of initial configurations. Curves with equal e and different JV 
collapse for e < e s . This shows that the dissipated work per monomer is an extensive function and AW ~ N 2 v. Again, 
this scaling behavior breaks down above the liquid-solid transition. 



dissipation without an actual reptation of the chain. We did not analyze the precise microscopic dissipation mechanism 
and therefore leave this issue for future work. 



Integrating the friction force, eq. (24 1, yields the scaling form of the dissipated work 



AW(v,e,N) 



F h .(v, e, N G )dx ~ i lG {e)a 2 N^ +1 v , 



(26) 



where we assume JVq ~ (L c — x)/(2a). In fig. 17 AW/N 2 is plotted versus v, the data with equal e < £ s and different 
JV collapse. This implies that the dissipated energy per monomer AW/N and also the internal friction is extensive, 
7 = 1, meaning that a finite fraction of the globule or the whole globule rearranges during pulling and not only a 
few monomers. This scaling behavior breaks down above the liquid-solid transition, which however does not imply 
that the friction becomes non-extensive but rather that linear viscous scaling is not valid anymore. It is also not clear 
whether the extensive scaling of the viscosity, defined here via the viscous force when pulling a chain segment out of 
the globule, should be expected to persist in the limit of very large globules, as the scaling suggests the viscosity in 
fact to diverge in this thermodynamic limit. More work on this is also needed. 
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4.3.3 Internal viscosity 



To extract the internal viscosity quantitatively, we fit the dissipated work to a linear form in v according to 

AW{v,e,N) = f(s,N) v , (27) 

and extract the prefactor r(e,N), which is apart from geometrical prefactors proportional to the friction coefficient 
r defined in eq. (24). As can be seen in fig. 14 AW ~ v for all velocities studied if e < e s . Consequently we include 
all velocity data in the linear fit. However, for s > e s a marked deviation from linear behavior is observed. There, we 
fit only to the slowest velocities, where the data still scales linearly with v. From the scaling r(e,N) ~ riG(e)N~' +1 
one sees that the relative internal viscosity follows as 



t?g(e) 



r(e,N) 
f(0,N) 



(28) 



which is shown in fig. [18] for annealed and un-annealed initial configurations. Note that in the limit of vanishing 
cohesion, e = 0, the globular viscosity t]q(e) is due to solvent friction effects only. The internal globular viscosity 
without solvent effects can therefore be defined by the difference r/Q(e) — t]q(0), but note that this definition is 
approximate since the solvent friction also depends weakly on the value of e since the globular structure changes with 
varying cohesive strength. By plotting the ratio of the viscosity in the presence and absence of cohesion, any residual 
numerical prefactors and the polymer length dependence in eq. ( 26 1 are eliminated and we are able to compare all 



different sets of parameters. The internal viscosities extracted from the non-equilibrium pulling sim ulations coincide 



for different N, which indicates that the scaling Ansatz r(e, N) ~ rjQ(e)N' y+1 works fine. In fig. 18 1 for the annealed 
pulling simulations we add data for the internal viscosity obtained from the equilibrium globule diffusion simulations, 
which are defined via ?7g( £ )/^g(0) = D /Dq(e) (the same data as already presented in fig. [6]d). The viscosity data 
from the equilibrium simulations are considerably lower for s < e s . This might in part be caused by an underestimate 



of the Rouse friction Dq defined in eq. (16), which enters the definition of the viscosity in the equilibrium simulations, 



or by additional dissipativc mechanisms in the non-equilibrium pulling simulations: sometimes it happens that the 
whole globule is moved through the solvent when an entanglement within the globule does not yield quickly enough, 
which causes additional solvent dissipation. Considering that the two ways of extracting the internal viscosity are very 
different, in terms of the geometry employed and the general setup (one being equilibrium, the other non-equilibrium), 
we consider the agreement sufficient at the present stage. 

In fig. |18| we compare our data to a simple model based on the effective friction coefficient of a single particle diffusing 
in the periodic potential U p — (9e/2) cos(ttx/ci) with amplitude Be. This potential mimics the energy landscape the 
chain monomers are experiencing as they are moving against each other during conformational reorganization processes. 
In writing the amplitude of the potential as 9e, we assume the corrugation strength to be proportional to the cohesive 
energy e. The numerical prefactor 9 is a scaling factor that may be viewed as a fitting parameter. The solution of this 
one-dimensional diffusion problem yields an effective viscosity |30j 



77(e) 

7,(0) 



J-n 



9e 
2k B T 



(29) 



3> 1 [62 . As can be seen from figs. 



in the stationary long-time limit. Io (z) is the zeroth order modified Bessel function with the asymptotic limits Io (z) 
1 + z 2 /4 for z < 1 and I (z) — 1 
from simulations, as defined in eq. (|28 1 



18 



the rescaled internal viscosity 
is reproduced quite well by the model prediction eq. ( |29[ ) with a fitting value 
9 = 1, which is shown by a solid line. Our results compare excellently to our previous results where smaller globules 
N < 100 have been considered |32) . 

In the solid regime for large cohesive energies e > e a , pronounced deviations between the equilibrium diffusional and 
the non-equilibrium pulling simulations appear, which are better appreciated when the whole data set is plotted on a 
logarithmic scale in fig. |18b . The pulling simulations give a much smaller viscosity when compared to the equilibrium 
diffusion simulations, which shows that the pulling simulations for e > e s do not reach the linear response regime, a 
fact that is independently suggested by the absence of scaling in fig. 14 



4.4 Summary for pulling on homopolymer globules 

The dissipated work per monomer scales linearly in the pulling velocity, indicative of viscous friction; the friction is 
also proportional to the number of monomers in the globule and scales extensive, thus a finite fraction of the globule 
contributes to the internal friction, fig. [17] We show that the dependence of the internal friction on the cohesive 
strength e is described well by the diffusion of a single particle in a corrugated potential. We extend our previous 
results [35] to significantly larger systems and are able to show that - below the liquid-solid transition - the history, 



i. e. the preparation of the initial structures, does not influence the scaling results, see fig. 18 
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Figure 18: Rescaled internal viscosity rjG(s)/r]Q(Q), eq. (28 1, as a function of the cohesive strength e for (a), (c) annealed 



(b), (d) un-annealed initial configurations. The viscosity ?7g( £ ) is obtained from linear fits to the dissipated work as 
a function of the velocity v. As the data points for different N coincide, all t(e, N) exhibit the same N dependence. 
The solid line is the prediction of the excess viscosity of a Brownian particle in a sinusoidal potential, see eq. ( 29 ) . 
In (a), (c) the viscosity ratio as obtained from our equilibrium simulations, see fig. [6]d, ??g( £ )/ 7 7g(0) = D / Dq(s) is 
shown in addition. 



5 Summary and conclusions 

Two different dynamical regimes exist for a globular homopolymer, a solid-like regime at large cohesive energies or 
low temperatures, characterized by a substantially increased force needed to unfold the globule and by pronounced 
fluctuations in the force extension traces, and a liquid-like regime at low cohesive energies or high temperatures. The 
critical cohesive energy e s at the transition between these two regimes depends on the size of the globule Nq , as shown 
in figs. |12| and |16| In the liquid regime the monomers inside the globule are rather mobile and the internal friction or 
viscosity in the globule can be extracted from simulations. The globule dynamics is studied by two different scenarios, 
(i) by considering the equilibrium diffusion of a globule relative to the linking straight chain sections, and (ii) by 
non-equilibrium stretching simulations. In both scenarios we find that the internal friction or viscosity is extensive, 
thus scaling linearly with Nq and increases with growing e until the liquid-solid transition is reached. The signature 
of the solid state is a vanishing globule diffusion on the simulation time scales, Dq = 0, for the equilibrium globule 
diffusion setup, while huge fluctuations in the force extension curves are seen for the non-equilibrium pulling setup. 
The solid state is characterized by very slow internal dynamics and no reliable estimates for the internal friction can 
be obtained from our simulations. 

Note that we define the internal viscosity in a rather specialized way, namely via the viscous force needed to pull a 
chain segment out of the globule (other definitions of internal viscosity are possible and useful in different contexts) . 
Our prediction that the friction of chain motion scales extensively with the globular size should be discussed in light 
of the classical reptation scenario [8]: Here, the friction coefficient (per monomer) is assumed to be independent of the 
globule size but rather the whole chain is assumed to move or reptate through the melt, also giving rise to an extensive 
scaling of the friction. This suggests that reptation is at the heart of the extensive scaling of the viscosity found in our 
simulations, but other mechanisms are also conceivable. Besides, it is not clear whether the extensive scaling of the 
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viscosity holds also in the thermodynamic limit and how a possible crossover is brought about. For a melt consisting 
of finite-length polymers, the chain size would constitute a possible crossover length, for the hypothetical limit of a 
single infinitely long polymer chain no such crossover length comes easily to mind. 

The most direct experimental realization of our system would be possible with a homopolymeric globule in a single- 
molecule setup. But there might also be implications for protein folding. Our results suggest that a quick collapse 
of a protein might result in a kinetically trapped and misfolded hydrophobic core, a so-called molten globular state, 
which would take a long time to fold into the native structure due to high internal friction. Since the friction scales 
extensively with the monomer number inside the globular core, the folding time becomes prohibitively long already 
for moderate globule size. From solvent-viscosity dependent measurements of folding times, the internal viscosity of a 
short a-helix forming peptide in the absence of solvent viscosity effects has been estimated to be of the same order as 



the solvent viscosity, i.e. t]g(£)/t]g(Q) — 2 in our notation [9]. The interaction energy we extract from fig. 18 follows 
to be of the order of e ~ 2.5kBT, not unrealistic for typical interaction parameters for residue-residue contacts. When 
speculating about the relevance of our results to protein folding, a few cautious remarks are in order: First, protein 
collapse is driven by the hydrophobic effect, which is quite involved due to the presence of water solvent, neglected 
in our treatment using an implicit solvent model with Lennard- Jones interactions. We would think that the scaling 
properties of the internal friction and internal viscosity will not change by adding explicit solvent, which however 
should be tested. Secondly, specific interactions between protein residues, sequence effects and the existence of a well 
defined native state will certainly alter the friction behavior as well. One lesson that might be learned from our work 
is that even for the relatively simple case of a homopolymeric globule, the concepts of internal friction and internal 
viscosity are far from trivial. 
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